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We study the interface tension of the 4-state Potts model in three dimensions using the Wang- 
Landau algorithm. The interface tension is given by the ratio of the partition function with 
a twisted boundary condition in one direction and periodic boundary conditions in all other 
directions over the partition function with periodic boundary conditions in all directions. With 
the Wang-Landau algorithm we can explicitly calculate both partition functions and obtain the 
result for all temperatures. We find solid numerical evidence for perfect wetting. Our algorithm 
is tested by calculating thermodynamic quantities at the phase transition point. 
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1. Introduction 

The interface free energy on a D-dimensional lattice, with size Li x L2 x • • • x Lo = Li x Vd_i, 
is obtained from the ratio of two partition functions 

Fs = -\og^-\ogU, (1.1) 

where Zp is a partition function with periodic boundary conditions and has periodic boundary 
conditions in D — 1 directions and a twisted boundary condition, which develop a topological ex- 
citation, in the direction labelled by 1. In other words, the topological excitation forced by the 
twisted boundary condition induces an interface in to the system, and the free energy of the inter- 
face is obtained by subtracting the free energy of the bulk (— logZp) out of it. The logLi term has 
been subtracted since the interface can form at any of the L\ points. The interface tension is then 

0= lim (1.2) 

L,^c» Vb-i 

The difficulty in calculating the interface tension arises from the fact that the partition functions 
are hard to calculate using traditional simulation methods, since they have exponential fluctuations 
in volume [1, 2]. Therefore, the standard approach has been to relate the interface tension to 
quantities that aie easier to determine with Monte Carlo simulations [3]. However, these methods 
generate large systematic and/or statistical en^ors. 

In [4] we have proposed a new method to calculate interface tensions. It is based on the Wang- 
Landau [5] algorithm. This algorithm enables us to calculate numerically the density of states of 
the system, and from this the partition function can be obtained as a numerical integral. As a test 
system we used the 3d 4-state Potts model. We were also able to provide solid numerical evidence 
for perfect wetting. 

2. Wang-Landau algorithm 

Consider a system with partition function 

Z = £exp(-i3£,-) = £^(£,)exp(-i3£,-)- (2-1) 

{'•} 

where {/} is a generic configuration and the density of states with energy £",-. The Wang- 
Landau algorithm [5] estimates directly the density of states g{Ei). The algorithm is based on a 
simple observation: If the state of a system is changed randomly, and the probability to visit a given 
energy level £ is 1 /g{E), a random walk in E generates a flat distribution across the energy levels. 
The algorithm consists of following steps 

1. Start with any lattice configuration and a nonzero density of states giE). 

2, Change a random site to a random value. 
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P = inin<;4^,U (2.2) 



3. Accept the change with a probabiUty 

= min i 

and if the change is rejected set Emv, = fold- 

4. Update g(£'new) and a histogram H{Eaew) 

) fg{Enew) , 
//(£„ew) ^ //(£new) + 1 ■ (2.3) 

5. Set Sold = £'new and repeat from 2 until g{E) has converged. Then reduce / — ^/f and set 
H{E) to zero. 

6. Repeat from 2 until / reaches some small enough value. 

The algorithm satisfies detailed balance only in the limit / — ?■ 1. In more detail, the transition 
probabilities satisfy 



^ P{E^E') 



1 



rME'^E), if g{E)>fg{E') 



' piE'^E), if giE') <g{E)<fg{E') . (2.4) 



\rME'^E), if g{E)<g{E' 



At / = 1 we obtain the detailed balance condition 

-^p{E ^ E') = -^ME' ^ E). (2.5) 
g[E) g{E') 

Hence, with / 7^ 1 the algorithm produces a g{E) that favors the large values over the small ones. 
However, in practice this does not seem to generate problems as long as the number of iterations is 
large enough, s.t., the statistical errors are larger than the systematic ones. 

As a convergence criterion we follow [6], and require that each energy level is visited 1 / y/lnf 
times since after this statistical enors are not expected to decrease. This criterion might suffer from 
systematic errors, which can lead the system to convergence to a wrong density of states [7, 8]. 
A safer criterion would be to require 1/ln/ visits to each energy level. However, for the level 
of accuracy requested in our simulations, the results obtained with these two criteria are identical 
within the statistical errors. See next section. 

3. Algorithm testing 

As a test case we use 4-states Potts model in 3d. The Hamiltonian of the model is 



('7) 



where / is the strength of the interaction, qi G [0,3], 5q.,q. is the &onecker delta function of the 
spin variables qi,qj on neighbor sites ij and the sum {ij) is over nearest neighbors. 
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Figure 1: Specific heat calculated using three different convergence requirements: flat histogram (left), 
1 / ln(/) visits required (middle) and 1 / -^/In/ visits required (right). The solid line is the result of the middle 
panel after 27 iterations. 
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Figure 2: Left: Normalized difference between the specific heat obtained with two different convergence 
criteria. The label a indicates the criterion requiring l/^/In/ visits at iteration 23 and b 1/log/ visits at 
iteration 27 . Right: Transition temperature calculated with three different methods. 



To obtain statistical error estimates for measured quantities we perform 20 separate runs witli 
twisted and untwisted boundary conditions. For each run we do 23 iterations, meaning that /min ~ 
1 + 1.3 X 10^^ during the last step (we start by setting / = 3). The simulations are done on lattice 
sizes L= 16,20,24,28,32,40,48, and 56. 

To test the convergence criteria and the systematic errors arising from the update, we measure 
the specific heat 




(3.2) 



We have used three different convergence criteria. The first requires fiat histogram, the second 
l/ln(/) visits and third l/y^ln(/) visits to each energy level. In Fig. 1 we have plotted the specific 
heat calculated with the different convergence criteria as a function of the number of iterations. At 
about iteration 20 all the three criteria provide convergence to the same value. In the left panel 
of Fig. 2, we have plotted the normalized difference of the specific heat obtained with the third 
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0.3143104(9) 




0.3143103(9) 




0.3143102(9) 


A: 


0.3143103(9) 




0.3143103(5) 


Pbbd 


0.3143041(17) 



Ae 


1.16454(16) 


AeMM 


1.16492(12) 


AeBBD 


1.16294(61) 



Table 1: Comparison of our results with previous studies for the transition temperature (left) and the latent 
heat (right). The result with label MM are from [10] and BBD from [9]. 



criterion at iteration 23 and the specific heat obtained with the second criterion at iteration 27. The 
differences are completely dominated by statistical eiTors, and hence we do not expect to have 
systematic en^ors due to the finite / or wrong convergence of g{E). 

To additionally test our algorithm we compute the critical temperature and latent heat and 
compare them with the results of [9] and [10]. Following [9], we use three different definitions 
for the critical j3 at finite volume. The first, {5^, is defined to be the value for which the canonical 
distribution P{E,j5) = g{E)e^^^ has two equal maxima, the second, j3^^, is the position of the 
central energy of the latent heat, and the third, is the position of the maximum of the specific 
heat Eq. (3.2). After exti'apolating to infinite volume these definitions agree with each other. See 
Fig. 2, right. 

The latent heat Ae can be obtained from the maxima of specific heat 

Cmax(L)=C+^(J8,)2(A^')2L^ (3.3) 

Our results and those of [9, 10] are collected in Table 1. The discrepancy with [10] is less than 
two standard deviations for all the quantities. The agreement is not as good with [9], but in this 
work smaller lattices are used. 

4. Interface tension 

We consider two different kinds of interface. The first, order-order interface, is between two 
regions of space that are in two different ordered vacua. Near the critical temperature, the interface 
tension assumes the form (see [11, 12]) 

Ooo(L) =000 + ^ + ^..., 

from which the asymptotic behavior can be extracted. In our simulations we have truncated the 
series to ^(L^^), since the corrections from higher order terms ai^e same the size of the statistical 
errors. For the fit we use only the points for which 
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Figure 3: Left: The order-order interface tension as a function of j8. Right: The probability distribution at 
the critical temperature. The maxima are normalized to one. 



The former condition is to reduce finite size effects and the latter one is dictated by the fact that 
Eq. (4.1) is only expected to hold at large distances. 

Near j3c, the behavior of Ooo can be parametrized as 

aoo(/3) = aoo(i3,)+a(/3-i3,)''- (4-2) 

This functional form represents the data well. A fit with 9 degrees of freedom has ;^^/9 = 0.11. 
The values are Ooo = 0.0249(6) and p = 0.76(4). See the right panel of Fig. 3. 

The other interface, the order-disorder interface, is an interface that separates an ordered state 
from the disordered state. It is defined only at the critical temperature for systems with a first order 
phase transition. The order-disorder interface can be calculated from the probability distribution. 
Namely, if P^ax is the peak of the histogram when the two maxima have equal height, and Pmin is 
the minimum between the peaks, the tension is given by (see [13, 14]) 

2aod = y^Iogf^) . (4.3) 

Our data for P are displayed in the right pane of Fig. 3. 

To obtain the asymptotic form we use a fitting function of the form. 

logL c? ca 

2cTod = - 2aod + ^ + ^ . (4.4) 

A fit of the date according to Eq. (4.4) yields 20do = 0.0252(4). The insertion of a term C3/L^ in 
the fitting function, justified by the use of the estimator (4.3), gives a compatible result. 
The ratio of order-order to disorder-order interface tensions gives the wetting angle 

;^ = cos0, (4.5) 
2aod 

where for perfect wetting = 0. See Fig. 4. Our result is a clear indication that perfect wetting 
holds for the three dimensional Potts model. This result has also been conjectured for the two- 
dimensional Potts model [15]. 
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